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ABSTRACT 

We present results from a multi-wavelength study of 29 sources (false detection prob- 
' abilities < 5%) from a survey of the Great Observatories Origins Deep Survey-North 

ly-j , field at 1.1mm using the AzTEC camera. Comparing with existing 850 |J.m SCUBA 

■ studies in the field, we examine differences in the source populations selected at the two 

wavelengths. The AzTEC observations uniformly cover the entire survey field to a 1-a 
depth of ~ ImJy. Searching deep 1.4 GHz VLA, and Spitzer 3-24 |i.m catalogues, we 
identify robust counterparts for 21 1.1mm sources, and tentative associations for the 
I remaining objects. The redshift distribution of AzTEC sources is inferred from avail- 

able spectroscopic and photometric redshifts. We find a median redshift of z = 2.7, 
somewhat higher than z = 2.0 for 850 /im-selected sources in the same field, and our 
lowest redshift identification lies at a spectroscopic redshift z = 1.1460. We measure 
jj] ■ the 850 |i.m to 1.1mm colour of our sources and do not find evidence for '850 |J.m 

' dropouts', which can be explained by the low-SNR of the observations. We also com- 

bine these observed colours with spectroscopic redshifts to derive the range of dust 
temperatures T, and dust emissivity indices (3 for the sample, concluding that existing 
estimates T ~ 30K and f3 ~ 1.75 are consistent with these new data. 

Key words: galaxies: high redshift - galaxies: starburst - galaxies: formation - 
infrared: galaxies - submillimetre. 
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1 INTRODUCTION 

Over the last decade observations at submillimetre and 
millimetre wavelengths (350-1200 vim) have been used 
to detect a population of luminou s (Lib. = 
1000 yim) > 10^^ Lq) galaxies (e.g. ISmail et al. 



Hughes et al] 1 19981: iBarger et all il998l:lE'ales et al.l 



Cowie et al.ri2002l : IScott et al.l |2002| : iBorvs et al.l 
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Serieant et al' '2003": 'Webb et all l2003l: IWans et al.' 


'2004; 


Grcvo ct al., 2004; 


Laurent et al. 2005: Coppin ct al. 


2005, 


2006; Knudsen et al. 2006: Bertoldi et al.l 200?!: iKhan et all 


2007.:.Scott et al.ll200&, Greve et al.l2008l:|Perera et al.ll2008l: 


Greve et al.l 


20081: Devlin et al.ll2009l. Austermann et al. sub- 



mitted). These objects, referred to as submillimetre galax- 
ies (SMGs), are thought to be high redshift analogues of 
the loca l Ultra Luminous Infra-r ed Galaxy (ULIRG) pop- 
ulation l|Sanders fc Mirabellll996l ). Their large luminosities 
and apparent lack of significant active galactic nuclei (AGN) 
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activity in most c ases (e.R. Hornschemeier et all 2000l: 
Bautz et al] I2OO0I: lAlmaini et al. l l2003l : [Alexander et al.l 
2OO5I : IPope et al.ll2008l l imply star-formation rates ^ 100- 



1000 Moyr-^ With orders of magnitude larger space den- 
sity at z > 1 than in the present-day Universe, it is presently 
believed that SMGs could represent an energetic early star- 
forming phase in the process that produces giant elliptical 
galaxies, and a significant fraction of th e total star-formation 
rate density at 2 ^ 2 fsee lBlain et al.li2002 . for a review). 

The identification of multi-wavelength counterparts to 
SMGs is hindered by the angular resolution of the current 
generation of submillimetre (submm) instruments (typically 
~ 10-20 arcsec), and the high surface density and faintness 
of counterparts in the optical/NIR, making unambiguous 
associations difficult. Significant progress has been made in 
the field by first searching for candidates in much lower sur- 
face density catalogues with higher astrometric precision, 
in particular using 1.4 GHz VLA interferometer maps, and 
deep 24 \im Spitzer observations. This method works both 
radio and mid-IR wavelengths, as in the submm, are biased 
toward the detection of star-forming galaxies: the radio syn- 
chrotron emission has a well-known correlation with the far- 
IR radiation that gets redshifted in to the observed submm 
band, and 24 yim samples primarily thermal emission from 
warmer dust in the vicinity of star-forming regions. With 
the much improved positional uncertainties of ~ 1 arcsec of- 
fered by these radio and mid-IR data sets, it is then possible 
to identify optical/ NIR counterparts provided that they are 
bright enough (e .g. Ivison et al.l 2002l : lchapman et al.ll2005l : 
iPope et al.ll2006l : llvison et al.ll2007i r 

In this paper we use this established procedure to iden- 
tify counterparts to SMGs detected in a 1.1 mm ma|j3 of the 
Gr eat Observatories O rigins Deep Survey-North (GOODS- 
N, H erera et all I2OO8D usi ng the Astro nomical Thermal 
Emission Camera fAzTEC. Iwilson et al.| [2008l. GOODS-N 
is one of several well-studied fields in the northern hemi- 
sphere that has the prerequisite radio data, as well as deep 
Spitzer coverage to identify counterparts. There is also an 
impressive collection of optical imaging {HST and ground- 
based), and optical spectroscopy for > 1500 targets with 
which to study the detailed properties of individual objects 
once their positions are known. 

Until recently, the most complete submm i mage towards 
GOODS-N was the SCUBA 850 urn map of iBorvs et all 
l|2003l ) (see also iPope et al.l I2OO5I : IWaU et al.lT2008l ) that 
was produced from a heterogeneous collection of data ob- 
tained by different groups with different observing modes 
jHughes et al] Il998l : [Barger et all I2OO0I : iBorvs et al.1 l2002l : 
ISerieant e t al.''2003l: IWang et al.ll2004h . This map produced 
a sample of nearly 40 sources, and was the subject of a de- 
tailed multi- wavelength study (|Borvs et al.l2004 : |Pope et al.l 
I2OO5I . [ 20061 ') . However, the spatially varying noise of the 
SCUBA map, combined with the desire to search for even 
higher-redshift sources that are expected to be more easily 
detected at longer wavelen gths due to the m ore favourable 
negative J f-correction ( e. g.lEales et al.ll2003l ), motivated the 
survey of IPerera et al.l (|2008h to uniformly map the en- 
tire area at 1.1mm. The GOODS-N AzTEC map covers 
245 arcmin'^ (matching the Spitzer coverage) , has an 18 arc- 



^ Map available at 'http: / /www. astro. umass. edu/AzTEC7] 



sec full-width half-maximum (FWHM) beam (compared 
with 15 arcsec for SCUBA at 850 ^.m) and reaches a uni- 
form RMS depth of 0.96-1.16 mjy beam"^ Note that there 
is also a map covering a sim ilar area made using MAMBO at 
1.2 mm l|Greve et al. I2OO8I ): those data have a smaller beam 
(11.1 arcsec FWHM), but slightly less uniform coverage with 

noise varying between 0.7-1.2 mjybeam"^. 

The 28 robust 1.1 mm sources identified in lPerera et al.l 

(|2008l ) were detected with significances > 3.8-a. In this pa- 
per we present potential counterparts for all of these sources, 
as well as one new object that was obtained by deblending 
the brightest peak in the map, AzGN 1, revealing a faint 
source that we label AzG N 1.2 (correspond ing to GN 20 and 
GN 20.2, respectively, in lPope et aLll2006l . Pope et. al sub- 
mitted). Note that the 1.1 mm deboosted flux densities given 
in this paper in Table IA3I have been corrected for Edding- 
ton bias and in many cases have signal-to-noise ratios (SNR) 
< 3-a: these values are the least biased estimates for the true 
flux densities, but do not reflect the robustness of the de- 
tections. The integrated negative tails of these distributions 
were used to estimate false detection probabilities, and a 
limit ^(5* < 0) < 5% corresponds to the 3.8-a threshold men- 
tioned above. Using extensive simulations the actual spuri- 
ous rate f or the entire sampl e was estimated to be 1 or 2 
sources m iPerera et all l|2008D . 

We find robust counterparts for 21 objects, which we 
define to be objects with false-identification probabilities 
P < 0.05 within 6 arcsec. We also provide tentative identifi- 
cations for the remaining sources, considering counterparts 
up to 10 arcsec away and 0.05 < P < 0.10 (Section[2]). These 
identifications enable us to report radio-IR spectral energy 
distributions (SEDs) using the VLA, SCUBA and Spitzer 
GOODS-N data. For the robust list, we identify spectro- 
scopic redshifts for 7 objects in the literature, and provide a 
combination of mid-IR and radio-(sub)mm photometric red- 
shifts for the remaining 13 sources (Section [3]). We compare 
our results with the existing SCUBA studies in this field to: 
(i) identify differences in the redshift distributions of sources 
selected at 850 p.m and 1.1mm; (ii) evaluate the effective- 
ness of searching for '850 ^.m dropouts' (objects detected at 
1.1mm but not at 850 n.m) as a means for finding higher- 
redshift SMGs; and (iii) probing the rest-frame distribution 
of dust properties of SMGs consistent with measurements in 
the two bands (Section 0}. 



2 COUNTERPART IDENTIFICATION 

2.1 Radio and mid-IR matching catalogues 

The radio and Spitzer catalogues th at we use to find coun- 
terparts are generally the same as inl Pope et al.1 l|2006l ). and 
we refer the reader to that paper for further details. The only 
significant update to their analysis is an improved 1.4 GHz 
VLA radio map, with a 50% reduction in the noise to ~ 4- 
5 fj,Jy RMS across the AzTEC cov erage region compared to 
that presented by iRichardsl (|2000|), and a bout 25 per cent 
deeper than the map used bv lPope et al.l 0006). The com- 
plete data set contains a total of 165.5 hr of VLA 1.4 GHz 
observations in A- (128.5 hr), B- (28 hr), C- (7hr), and D- 
(2hr) configuration. These data were combined, reduced, 
and imaged using AIPS. Full details of this analysis will be 
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Figure 1. The relative coverage of data sets in GOODS-N used in this paper. The greyscale indicates the RMS noise in the AzTEC 
1.1mm map. The soUd black contour corresponds to a noise of 1.16 mjy in this map, and is the region within which AzTEC sources 
(white numbered triangles) were extracted. The white dashed contour indicates the SCUBA 850 |am coverage with a noise less than 
10 mJy. The black dashed contour shows the MIPS 24 \xm coverage, and the dot-dashed lines the IRAC 3.6 \xm coverage. Finally, the thin 

dotted contours indicate the surface density of VLA 1.4 GHz sources from the 4-ct catalogue (0.5, 1.0, 1.5 and 2.0arcmin^^) measured in 
apertures with a radius of 0.1 dcg. For reference, the half-power radius of the radio map noise is about 0.25 deg. The holes seen in these 
contours centred over 12:34:52 +62:03:41, 12:35:38 +62:19:32, and 12:38:48 +62:23:09 are regions of the map that were excised due to 
sidelobe structure caused by especially bright sources. 



presented in Morrison et al. (in preparation). While the fi- 
nal product of that paper will be a S-cr catalogue, we have 
produced two deeper catalogues for use in this work: (i) 
an approximately 4-cr catalogue with a surface density of 
1.80 arcmin"^; and (ii) a fainter 3-cr radio catalogue with 
a surface density of 3.52 arcmin"^ that initially contains a 
much larger fraction of spurious sources from which we se- 
lect only objects that are coincident within 1 arcsec of 3-cr 
detections in the 24 [im catalogue. In both cases a radio 
catalogue of the given significance is constructed using the 
AIPS source extraction task SAD in signal-to-noise mode 



(using the uncorrected peak flux densities) which uses an 
RMS map estimated from the task RMSD. After SAD de- 
tects the sources over a region encompassing the AzTEC 
data, it then fits Gaussian components to these sources, 
and in the process applies corrections for both the band- 
width smearing and primary beam attenuation. The final 
catalogues report these corrected flux densities, including 
any spatially resolved structure that increases the extent of 
the fitted Gaussians, without any additional 'by-hand' re- 
moval of problematic sources. The catalogues are therefore 
expected to contain a number of false positives, and begin 



4 Edward L. Chapin et al. 



1.0 




101 



w 

'c 10° 

6 

O 

^ 10-1 

m 

A 

10-2 



10 100 1000 

s. [^Jy] 

Figure 2. Integral source counts in tlie radio and mid-IR match- 
ing catalogues. The primary catalogue is the i-a radio catalogue 
and its integral source counts within the solid black contour in 
Fig.[T]are shown by a solid line. A deeper 3-(t radio catalogue was 
produced by the intersection of a raw 3-(t radio catalogue with the 
3-(T 24 |am catalogue. The counts from the two raw source cata- 
logues are shown by the dotted and dot-dashed lines respectively, 
measured within the area common to the solid black and black 
dashed contours from Fig. [T] While the raw S-cr radio catalogue 
has a huge increase of sources at Si. 4 < 25 ^tJy, indicating many 
spurious detections, the requirement of a 24 |xm counterpart re- 
jects many of these sources, such that the remaining counts at 
fainter flux densities follow the trend established at the bright 
end, but extending to a fainter limit than the 4-(T catalogue. 



to suffer incompleteness, at their respective flux density lim- 
its. The 3-cr catalogue covers a slightly smaller area, as it is 
limited to the 24 \im coverage. For several AzTEC sources 
around the edge of the map we therefore use only the A-g 
radio catalogue to make identifications) ; see Fig [T] for the 
relative coverage of each data set. 

In Fig. [2] the integral source counts for the radio and 
mid-IR catalogues are shown. The horizontal axis has been 
extended to sufficiently faint flux densities to show the point 
at which each catalogue becomes incomplete (flat integral 
counts). Clearly the raw 3-cr radio catalogue contains many 
spurious detections at 51.4 25 /xjy, as the counts diverge 
steeply from the trend at brighter flux densities traced by 
both the 3-(T and A-a catalogues. However, the requirement 
of a 24 \im counterpart for each radio source drastically re- 
duces the number of candidate detections in this fiux density 
regime, effectively extending the faint counts to ~ 15-20 /iJy 
from the ~ 20-25 /.iJy achieved in the 4-(t catalogue. We 
note that while this technique enables us to extract fainter 
sources from the radio map (rejecting many of the spuri- 
ous noise peaks) , the requirement of 24 fi.m emission could 
lead to incompleteness (i.e. real radio sources that are not 
detected with MIPS), as the 24 (i.m channel samples rest- 
frame spectral features (e.g. polycyclic aromatic hydrocar- 
bon emission and silicate absorption) that are not directly 
related to the radio emission. Nevertheless, it will be shown 
in the following sections that even with these potentially in- 
complete radio catalogues we identify a significant fraction 
of the 1.1 mm sources. 

For simplicity, a fixed search radius was used to identify 
potential counterparts in the two catalogues, rather than 
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Figure 3. The fraction of AzTEC sources with robust 4-cr 
1.4 GHz (P<0.05) candidate identifications as a function of 
search radius (solid black line) compared with the fraction of 
sources with any candidate identifications (no cut on P, solid 
grey line). For reference the black and grey dashed lines indicate 
the spurious counterpart detection rates using random positions 
with and without the P < 0.05 cut, respectively. A search radius 
of 6 arcsec (vertical dotted line) is used since it gives the peak 
robust identification fraction. It also happens to be where the 
random fraction corresponds to 5 per cent. We chose this 6 arcsec 
radius to search for all candidates, even in the 3-(t radio cata- 
logue. In the event that no robust counterparts with P < 0.05 
can be identified within 6 arcsec, the seaxch radius is extended 
out to 10 arcsec to find tentative identifications. 



a variable radius as a function of AzTEC SNR (see e.g. 
Ilvison et al.ll2007l ). There is also a relatively small dynamic 
range in AzTEC fiux densities, and with fairly constant 
noise, so that most of the positional uncertainties would be 
similar in any case. Given a search radius and list of po- 
tential counterparts, the probability that a given candidate 
is a rando m association, 'P', is c alculated following the pre- 
scription of lDownes et al.l (Il986l ). a method that is now used 
almost ubiquitously in the submm literature. This technique 
accounts for the surface density of sources in the matching 
catalogue as a function of brightness: if two potential iden- 
tifications of different brightnesses lie at the same distance 
from the AzTEC source the rarer brighter object will be as- 
signed a lower value of P. We have used the raw measured 
integral source counts in the matching catalogues (Fig.[2ll for 
this calculation, rather than a model. Since the radio cat- 
alogues undoubtedly contain a number of spurious sources 
near the detection thresholds, this procedure will naturally 
account for them by down-weighting the robustness of pro- 
posed faint identifications where the surface density is much 
greater. 

We also investigated the use of spatially-varying mea- 
surements of the counts. Since the noise depth of the radio 
map falls off towards the edge of the AzTEC coverage, so 
do the counts in the radio catalogues. We chose an aperture 
with radius 0.1° (significantly smaller than the ~ 0.25° half- 
power radius of the VLA map noise) in which we measured 
the integral radio counts centered over each position. For the 
4-0" radio catalogue this resulted in a variation of the total 
surface density ranging from about 2.4arcmin~^ at the cen- 
tre of the map, to 1.0-1.5 arcmin"^ along the edges (Fig.[T]). 
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We experimented with measurements of P using these mod- 
ified counts at the locations of each AzTEC source. With the 
increased source density toward the centre of the map, val- 
ues of P are increased slightly, and conversely values of P are 
decreased slightly toward the edges. However, the changes in 
individual values of P are generally < 20%, and the only ef- 
fects of using this calculation on the final list of robust identi- 
fications would be to add AzGN 2, as P would drop to 0.050 
from 0.051 for the single radio source within 6arcsec, and 
similarly, the second radio source near AzGN 18 would be 
added, with P dropping to 0.047 from 0.051 (thi s object was 
already identified as a radio double, GN 38, in IPope et al.l 
l2006l . and the photometric redshift estimates are consistent 
with them being at the same distance) . Since the differences 
are small, and the simpler calculation gives a slightly more 
conservative list of potential counterparts, we elected to use 
only the radio and mid-IR number counts averaged across 
the entire region of AzTEC coverage presented in Fig. [2] Av- 
erage source counts were also used by I Pope et al.l (|2006 ) and 
iGreve et all l|2008l ') in their calculations of P to find SMG 
counterparts in this field. 

2.2 Choice of search radius 

A search radius was chosen to provide a reasonable level of 
completeness, while minimizing the number of false identi- 
fications. In Fig. [3] we plot the fraction of AzTEC sources 
with at least one such counterpart in the 4-ct radio catalogue 
as a function of search radius with P < 0.05 (solid black 
line). The function initially rises as most true counterparts 
are eventually detected, reaches a peak of 62 per cent at 
6-7 arcsec, and then drops as potential identifications again 
become improbable due to the increasing chance of a spu- 
rious detection with such a large search radius. In contrast, 
if all radio sources within the search radius are considered, 
the total fraction with potential counterparts continues to 
grow (solid grey line). To demonstrate the effectiveness of 
cuts on P to lower the false identification rate, we repeat 
these calculations using random positions. Using the cut on 
P (dashed black line) results in a plateau at the expected 
fraction of ~ 0.05, coincidentally, also at a search radius of 
6 arcsec. Without the cut, the spurious identification frac- 
tion (grey dashed line) continues to rise to ~ 100 per cent 
by a radius of 60 arcsec. 

The chosen search radius of 6 arcsec is smaller than 
those used for SCUBA sources, typically in t he range 7- 



8 arcsec for objects w ith similar SNR (e.g. Ilvison et ahl 
2002;J\^febb ct al. 2003,; iBorvs et al.ll2004l : IPope et al.ll2006l : 
Ivison et al.ll2007l). despit e AzTEC having a slightly larger 



beam. Borvs et al.l (|2004l ) used a different method to esti- 
mate a search radius for counterparts to SCUBA sources 
in GOODS-N, finding a value of 7 arcsec to be appropri- 
ate. We repeated their analysis with our data, however, and 
concluded again that we should use 6 arcsec. This test sug- 
gests that the smaller resulting search radius is a property 
of the data, rather than the method we used to calculate 
it. As a further consistency check, recent SMA follow-up 
of AzTEC sources detected in other fields using the same 
pointing model as that e mployed here has sho wn excellent 
agreement. For example, lYounger et al.l ()2007l ') found posi- 
tional uncertainties < 4 arcsec between SMA and AzTEC 
centroids of millimetre sources in the COSMOS field de- 



tected with generally higher SNR than the objects discussed 
in this paper. 

Finally, we note in Fig. |3] that the fraction of AzTEC 
sources with any 4-cr radio candidates beyond 6 arcsec (no 
cut on P, grey solid line) continues to grow significantly 
faster (from ~ 65% to 80%) than the spurious rate (grey 
dashed line, from ~ 5% to 10%) to a search radius of 8 arc- 
sec. This suggests that ~ 70% of the entire AzTEC sample 
has counterparts in the A-a radio catalogue, of which 10%, 
about 3 objects, are within 6-8 arcsec of the AzTEC cen- 
troids. There is a radio source within 47 arcsec of every 
AzTEC centroid, although the additional objects encoun- 
tered beyond ~ 10 arcsec are almost certainly chance align- 
ments. 



2.3 Radio and mid-IR identifications 

We first search for counterparts in the 4-a radio catalogue 
with P < 0.05, and if none are found within 6 arcsec we 
proceed to search in the 3-cr catalogue using the same radius 
and cut on P. This procedure gives us a robust sample that 
will subsequently be analyzed in detail. In cases where no 
such counterparts can be found, we relax the search to iden- 
tify more tentative counterparts out to a radius of 10 arc- 
sec, and/or P < 0.10. Since these search parameters are 
expected to result in a much higher fraction of chance align- 
ments, this extended catalogue is not used to measure prop- 
erties of the general 1.1mm galaxy population, and is sim- 
ply included for completeness (although statistically speak- 
ing most of these tentative identifications are probably cor- 
rect). In one case {AzGN 27) neither radio catalogue yields 
a source within 6 arcsec, nor are there any counterparts with 
P < 0.05 out to 10 arcsec, so we search in the MIPS 24 yun 
catalogue (using the counts in Fig. [2] to calculate P), find- 
ing two potential identifications (however, we are still able 
to estimate radio flux densities by performing photometry in 
the VLA map at the MIPS positions). While possibly con- 
taining the correct counterparts for this particular object, 
the 24 (j.m catalogue is not, in general, as useful as the ra- 
dio catalogues, due to a significantly greater surface density 
and hence chance of identifying random interlopers (see also 
Ilvison et aLll2007l ). 

In total this procedure yields at least one counterpart 
within 6 arcsec for 26/29 AzTEC sources, 22 of which are 
robust with P < 0.05. Of the robust identifications, 18 were 
found in the 4-cr radio catalogue, 3 in the 3-a radio cata- 
logue, and 1 in the MIPS 24 \im catalogue. Of the tentative 
associations within 6 arcsec, 1 was identified in the 4-a radio 
catalogue, and 3 in the 3-a radio catalogue. Finally, all of 
the remaining 3/29 AzTEC sources have at least one ten- 
tative counterpart in the range 6 < r < 10 arcsec in the 
4-0" radio catalogue. These results are summarized in Fig. |4] 
and Table IX2I The spectral energy distributions for the pro- 
posed counterparts are given in Table IXSl Additional notes 
for each source can be found in Appendix [B] 



2.4 False identification rate 

We can estimate the number of spurious identifications for 
the 26 potential counterparts found within 6 arcsec by sum- 
ming their P values, which gives 0.85. This implies that 
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Figure 4. 30 arcsecxSO arcsec postage stamps showing AzTEC counterpart identifications in GOODS-N. The left panels show the 
VLA 1.4 GHz map stretched between — 20 /ijy (black) and +30 /xjy (white). The contours indicate fractions 0.1, 0.5 and 0.9 of the peak 
SNR in the AzTEC (white) and SCUBA (black; not always available) postage stamps. In addition, AzGN 1, AzGN 1.2, AzGN 14, 
and AzGN 22 have had the effects of nearby blended sources removed. The un-deblended AzTEC map contours are shown with dashed 
white lines for comparison. Cross-hairs indicate the locations of potential counterparts from the catalogues in which they were originally 
identified: white for the 4-(T radio catalogue; yellow for the 3-(T radio catalogue; and green for the MIPS 24 [jm catalogue. Long cross-hairs 
correspond to IDs with P ^ 0.05, and small cross-hairs for P > 0.05. The right panels are false-colour images constructed from IRAC 
3. 6 (blue), 4.5 (green ) and 5.8 ^mi (red) exposures. Squares indicate the locations of radio/IR counterparts to SCUBA sources proposed 
in lPope et al.l |2006|). The solid and dashed circles indicate search radii of 6 arcsec and 10 arcsec, respectively. Approximately 1 or 2 of 
the AzTEC sources are expected to be false positives. 
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of those 26 AzTEC sources, about one of the identifica- 
tions within 6 arcsec is expected to be spurious (noting that 
several sources have multiple proposed identifications). To 
understand how to interpret this result in terms of overall 
completeness, we consider several factors. First, the AzTEC 
source list is expe cted to have ~ 1-2 spurious detections 
(|Perera et al.ll200^ . Second, due to positional uncertainties, 
some of the true counterparts will lie beyond 6 arc s ec. W e 
adopt the radial offset distribution of llvison et al. I (|200'if ). 
rexp(-r^/2cr^), with cr ~ 0.6 x FWHM/SNR, which as- 
sumes a symmetric Gaussian beam and uncorrelated map 
noise. The cumulative distribution of this analytic PDF re- 
su lts in a shape very similar to the numerical simulations 
of IScott et aiTl|2008fi for AzTEC sources in the COSMOS 
field. Taking FWHM = 18 arcsec, and the SNR for raw map 
fiux densities (before deboosting), we would only expect to 
encounter counterparts within 6 arcsec for 27.5/29 sources 
on average if they were all real (neglecting positional un- 
certainties in the matching catalogue). However, we would 
confidently expect to find all of the objects within 10 arc- 
sec. Since we do not know which (if any) of the AzTEC 
sources are false positives, we simply apply this fraction to 
the expected number of real sources calculated above, and 
find that there should be ~ 26-27 real sources with true 
positions within 6 arcsec of their 1.1mm centroids. This ex- 
pectation is consistent with our identification rate within 
6 arcsec of ~ 25-26 out of 29 sources, although we stress that 
this statistical argument does not necessarily imply that the 
unmatched sources are spurious. We find P values less than 
0.05 for only 21 of the sources encountered within 6 arcsec 
(excluding AzGN 14 as noted below). While it may be the 
case that most of the remaining 5 sources are in fact asso- 
ciated with the 1.1mm objects, it is also possible that the 
true counterparts are simply fainter in the radio than the 
catalogue limit. We cannot distinguish between these two 
cases in the present study given the positional uncertainties 
in the AzTEC centroids. 

As a final warning, as with all studies that use P to 
evaluate chance alignment probabilities, there is an under- 
lying assumption that the matching catalogues are spatially 
unclustered. Two examples of ways in which this condition 
could be broken are an increased surface density of sources 
in the vicinity of SMGs due to multiple catalogue entries be- 
ing associated with the same physical structure (such as a 
galaxy cluster), or foreground lensing of background objects 
(an effect which is in fact c ommonly used to identify faint 
SMGs, e.g. lSmail et al.ll2002l ). In both of these cases P would 
be biased low. We do not attempt to correct for these effects 
in this work, but we alert the reader that evidence for such 
cases in GOODS-N will be discussed in the following sec- 
tions: the particularly complicated identificat ion of a coun- 
terpa rt for AzGN 14 (also known HDF 850.1. lHughes et al.l 
1 19981 ) which caused us to drop it from the analysis in this pa- 
per; and the potential presence of a protocluster at redshift 
z ~ 4. 



2.5 Comparison with SCUBA identifications 

Since 12 of the 29 objects discussed in this paper were also 
detected by SCUBA (see Table E2)) and i dentified in the ra- 
dio and mid-IR using similar techniques l|Pope et al.ll2006l : 
IWall et al ] |2008l ), it is useful to compare proposed identifica- 



tions to see how the new AzTEC positions and deeper radio 
catalogues affect the results. We exclude AzGN 14/GN 14 
{HDF 850.1) from this comparison (and most of the re- 
maining analysis in this paper) as its true counterpart has 
been under debate for some time due t o the suspected ob- 
scuration by a for eground elliptical fsee lDunlop et al.]|2004l : 
ICowie et al.l |2009| . and notes in Appendix [Bj . Of the re- 
maining 11 overlapping sources, we propose identical coun- 
terparts for 8 of those objects. For AzGN 1.2 we find 
that both the proposed counterpart of IPope et all (|2006l ) 
{P — 0.005 in this paper) and a se cond fainter radio ob- 
ject (with P = 0.037, also noted bv lDaddi et al.ll2008t ) are 
both robust identifications by our definition. In another sim- 
ilar case, only one object from a radio double identified by 
iPope et all l|2006l ') for AzGN 18 is strictly a robust coun- 
terpart (P = 0.032 in this paper), while the second radio 
source misses the cut, with P = 0.051. The only object 
for which we propose a completely different counterpart is 
AzGN 11 /GN 27, which was classified as 'tent ative' in the 
SCUB A map: the ACS/IRAC identification from lPope et all 
(200C|) is absent in the 1.4 GHz map, and we instead propose 
a radio source that lies slightly to the north, with P = 0.027. 



3 REDSHIFT DISTRIBUTION 

Some recent surveys at 1.1-1.2 mm claim to detect 
higher redshifts th a n SCUBA survey s at 850 ^im (e.g. 
lYounger etHI l2007l : iGreve et al1l2008 h. while others find 
redshift distributions that are indistinguis hable, possi- 
bly due to small samples sizes (e.g. iGreve et all |2004| : 
iBertoldi et"aLll2007l ). For this AzTEC survey we more accu- 
rately quantify any differences using greatly improved red- 
shift information, and comparing directly to the SCUBA 
results in this field using the same methodology. While the 
uncertainty in the 1.1 mm distribution derived from our data 
is large, due to the relatively small sample size, and cosmic 
variance resulting from the area of GOODS-N, this differ- 
ential measurement yields a useful comparison between the 
two bands. 

A number of groups have obtained spectroscopic red- 



shifts in GOODS -N (e 
|2004l: IWirth et al.l 



2006; Bar ger et al. 



Cohen et al? '2000'; 'Cowie et al" 



2004 : Chapman et a l. 2005; Rcd dv et al 
200a : (Pope et al.ll200aiDaddi et al.ll2008l . 



,2009 , Stern et al. in preparation). We found spectroscopic 
redshifts for 10 of our proposed AzTEC counterparts in 
these publicly available data-sets (see Table 1X2)) . However, 
one of those redshifts {AzGN 8) corresponds to the least 
favourable counterpart within the search radius (see the 
discussion for this source in Appendix [Bj . Another simi- 
lar case is AzGN 27 for which a spectroscopic redshift has 
been obtained only for the more distant of two potential 
counterparts. Finally, two radio sources that appear to be 
associated with the single object AzGN 7 lie at redshifts 
z — 1.996 and z = 1.992, and we assign a single redshift of 
z = 1.994, which is sufficiently precise for the purposes of 
this paper. Therefore our sample of 21 sources with unam- 
biguous identifications contains only 7/21 sources with spec- 
troscopic redshifts. This fraction is considerably lower than 
the 15/2 spectrosc o pic r edshifts for robust counterparts 
from the lPope et al.l l|2006l ) SCUBA sample (also excluding 
HDF 850.1), including the two new redshifts for GN 20 and 
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GN 20.2 from iDaddi et all (|2008l ). the redshift for GN 10 
from lDaddi et aL 1I2OOE ) and two additional Spitzer IRS red- 
shift sfownlPoE^ et all j2008l ). However, we are not surprised 



at this lower rate since we rely on archival data for the red- 
shifts of counterparts to new AzTEC sources, whereas many 
of the spectroscopic redshifts for SCUBA sources were ob- 
tained using targeted follow-up of proposed identifications. 

For sources without a spectroscopic redshift we first 
searched for optical photometric redshift estimates. As none 
were found, we instead employed two photometric redshift 
calculations using longer wavelength data. The first, Zir, is a 
simple function of the Spitzer photometry with coefficients 
derived from fits to SCUBA sources i n GOODS-N with spec- 
troscopic redshifts (|Pope et al.ll2006h . Although this method 
does not assume any particular SED, it benefits from the 
1.6 [im stellar bump that produces a strong characteristic 
feature in the observe d IRAC 3.6-8.0 \im bands for sources 
at re dshifts 1 z <^ 4 l|Simpson fc Eisenhardtll 19991 : ISawickil 
|2002| ). Such an empirical calculation may provide less bi- 
ased results than fitting spectral templates to the data since 
there are degeneracies between the derived redshift and as- 
sumptions about the starburst pr oducing the stellar bump 
(see discussion in lYun et all 120081 ) . While the residuals for 
this functional fit are relatively small (with a maximum 
Az = 0.4), no uncertainties are provided in iPope et al.l 
(|200^ ) for the remaining sources. However, a similar pho- 
tometric redshift estimator was derived by iWilson et al.. 
l|2008D . and a comparison with spectroscopic redshifts for 
SMGs from several fields (including SCUBA sources from 
both GOODS-N and SHADES) finds that a 1-cr error enve- 
lope Az = 0.15(1 -|- 2;) is a reasonable uncertainty estimate 
for 15 SMGs at reds hifts ^ z <, 3. We ha ve compared 
the lPope et"ai] l|2006l ') and I Wilson et al.1 (|2008| ) photometric 
redshift formulae for our data and find that of the 7 ro- 
bust identifications with spectroscopic redshifts both meth- 
ods provide estimates consistent with the spect roscopic mea- 
surem ents for the 3 sources at 2; < 3, within the lWilson et al.l 
l|2008l ) error envelope. However, both estimate s are b iased 
low at z > 3, more so using the IWilson et al. I l|2008l ) red- 
shift estimator. We also checked the scatter between the two 
methods for all of the robust identificati ons finding th a t they 
both gave answers compatible with the I Wilson et al] (|2008l ) 
uncertainty estimate. This bias and scatter are unsurprising 
as both formulae were fit to SMGs with spectroscopic red- 
shifts 2: ^ 3. In this work we assume the 1-a uncertainties 
are also Azii — 0.15(1 -I- 2), but warn the reader that the 
redshifts of more distant objects are probably systematically 
underestimated with this technique. 

The second photometric redshift indicator, 2rm, uses the 
radio and (sub)mm fiux densities fit to templates of local 
galaxies as suming the radio-IR correlation holds at hi gh red- 
shift (e.g. ICarilli fc Yunlll999l : lAretxaga et al.ll2007l ). This 
method provides the only redshift estimates for a handful of 
sources around the edges of the AzTEC map where there is 
no Spitzer or optical coverage from the GOODS survey (the 
entire AzTEC survey area overlaps with the 1.4 GHz data). 
O ur redshifts are calculated using the same methodology 
as lAretxaga et all (200 7) and summarized in Table 1X21 We 
note that the quoted 6 8 per cent confidence intervals are 
theoretical estimations; lAretxaga et al.l (|2007h checked the 
scatter between photometric and spectroscopic redshifts for 
a sample of SMGs with radio and submm data of similar 
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Figure 5. The redshift distribution of the 18 robustly identified 
1.1mm sources in GOODS-N (black histogram) with redshifts, 
adopting unique identifications from Table IA2I with the small- 
est P values for each source. The grey histogram is the redshift 
distribution of 20 similarly robust SCUBA 850 [am sources from 
IPope et al.l l|2006h , exclu ding GN 14, and updating some spectro - 
sco pic redshifts based on lDaddi et al.l 1I2OO8I) . iDaddi et al.l teOOOl ') 
and lPope et al.l l|2008l ) . There is not an integral number of sources 
in each bin because the uncertainties in photometric redshifts 
have been included. The results of K-S and M-W tests indicate 
chance probabilities pxs = 0.05 and Pmw = 0.01, that both his- 
tograms were drawn from the same parent distribution (although 
they can be larger when individual redshift uncertainties are in- 
cluded, see Section |3}. For comparison, the dashed red and green 
dotted histograms show the distribution of spectroscopic redshifts 
for AzTEC and SCUBA sources respectively. The excess number 
of objects in the 4.0 < z < 4.5 bin may be m embers of a high- 
redshift protocluster jPaddi et al.ll2008l. [ioool) . 



quality to GOODS-N finding an empirical symmetric 1-a 
scatter of Azrm = 0.8. We also adopt this uncertainty for 
consistency with the empirically measured uncertainties for 

Ziv. 

In Fig. [5] we show the AzTEC redshift distribution of 
the 18/21 robust [P < 0.05) identifications (black line) 
for which there are redshifts (spectroscopic or photomet- 
ric). The remaining 3 objects (^2;GA'' 15, AzGN 19, and 
AzGN 25) have only lower limits on their redshifts, and have 
been excluded from the histogram and subsequent analysis. 
For the objects lacking spectroscopic redshifts we use the 
photometric redshift estimates, including their uncertainty 
distributions to divide them amongst several bins if neces- 
sary. When both 2ir and 2rm are available, we calculate the 
variance-weighted mean redshifts. We feel we are justified in 
doing this because the two estimators depend on data and 
spectral features at completely different wavelengths, and 
are therefore independent of each other. We then calculate 
the distrib ution for 20 simila rly robust 850 (im sources (grey 
line) from lPope et al.l l|2006f ) in the same way for compari- 
son. Note that, as in this work, we have excluded HDF 850.1 
{GN 141 AzGN 14) from then redshift distribution. In ad- 
dition we have update d redshifts for several sources using 
new spectroscopic data l|Pope et al ] |2008l : lDaddr et al. I l2008l . 
I2OO9I) . 

The median redshift of the AzTEC sample is 2 = 2.7, 
with an interquartile range of 2.1-3.4. In contrast, the me- 
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dian of the SCUBA sample shown here is z = 2.0 with an 
interquartile range 1.3-2.6. For reference, the 14 SCUBA 
sources with spectroscopic redshifts (green dotted histogram 
in Fig. [5]) has a median z — 2.0, and the 7 AzTEC sources 
(red dashed histogram) z — 3.19, both in good agreement 
with the full distributions despite the small numbers of ob- 
jects. It is also worth noting that the spike of spectroscopic 
redshifts in the 4.0 < z < 4.5 bin seen at both wavelengths 
was found via ta rgete d follow-up t o SC UBA sources by 
iDaddi et all (|2008l ) andjDaddTelaD (l2009l ) that are thought 
to be members of a proposed protocluster at z ~ 4. Al- 
though not i nclud ed in the distribution, a recent study by 
ICowie et al.] l|2009l ) suggests that HDF 850.1 may also be a 
member of this high-redshift structure. 

We note that while the AzTEC sample appears to lie 
at slightly higher redshift than the SCUBA sample, we have 
had to rely more heavily on h ighly uncertain photometric 
estimates than in iPope et"aL I (l2006l) . However, the bias is 
likely to be toward lower rather than higher redshifts due to 
the nature of Zir. 

Next, taking the SCUBA and AzTEC redshift distribu- 
tions at face value (including photometric redshifts) we have 
used the Kolmogorov-Smirnov (K-S) and Mann- Whitney U 
(M-W) non-parametric tests for assessing how different they 
are. These tests are fair since both samples were drawn from 
the same region of space, and no extra uncertainty needs to 
be included to account for cosmic variance. Both methods 
operate on discrete samples so we first assign the mean red- 
shift to each object from its uncertainty distribution. The 
K-S test, which is sensitive to more general differences in 
the distributions (both the central values and tails) , gives a 
chance probability pKS =0.05 that both samples were drawn 
from the same parent redshift distribution. The M-W test, 
which is mostly sensitive to differences in the central val- 
ues of the distributions, gives a smaller chance probability 
Pmw = 0.01. However, we note that the uncertainties for ob- 
jects with photometric redshifts can be as large as Az ~ 0.8, 
comparable to the width of the entire population. To evalu- 
ate the spread in K-S and M-W probabilities that are con- 
sistent with our sample, we generate 10,000 mock samples 
at each wavelength drawing individual redshifts at random 
from the uncertainty distributions for each object. We find 
that 68 per cent of the time we obtain values pKS < 0.15 
and Pmw < 0.04. These tests show that, even with large in- 
dividual uncertainties, the shift to higher redshifts at 1.1 mm 
compared to 850 \im appears to be statistically significant. 

The SCUBA sample consists of a broader dynamic 
range in flux density than the AzTEC sample, due to the 
varying map depths, and the fact that different redshift pop- 
ulations may be present in t he deep and shal low regions of 
the map (|Pope et all l2006l : IWall et all 120081 '). It therefore 
may be the case that at least some of the differences be- 
tween these distributions are a result of a depth rather than 
wavelength selection effect. 

One potential concern with this comparison is that, due 
to a bias to higher r edshifts at fainter flux densities (e.g. 
IChapman et al.ll2005l ). the deeper radio catalogues used for 
matching in this survey simply detect mo re distant poten- 
tial counterparts than in lPope et al.l (|2006l ). We checked the 
distribution of radio brightness with redshift for our sample 
and found that the 6 faintest proposed radio counterparts 
lie in the redshift interval 2 < 2 < 3. Removing them, while 



broadening the remaining redshift distribution slightly, does 
not shift the median appreciably. However, if we remove 
sources with even brighter radio flux densities, we in fact be- 
gin to bias the sample to /lig/ier redshifts. Combined with the 
fact that we find most of the same counterparts for sources 
that appear in both the AzTEC and SCUBA surveys (Sec- 
tion [23]), we conclude that the intrinsic rest-frame scatter 
of radio luminosities in SMGs dominates any differences in 
the radio properties of 850 (xm and 1.1 mm selected samples. 

The lowest redshift that we find is AzGN 23 at z = 
1.146. This demonstrates the ability of mm- wavelength sur- 
veys to effectively select galaxies at z > 1, with little con- 
tamination from nearby objects. Assuming that our iden- 
tification procedure and redshift estimates are correct, and 
given the completeness of our survey, there is therefore lit- 
tle room for a significant tail to extremely high redshifts. 
Since the negative A'-correction at 1.1mm could in princi- 
ple enable us to d etect SMGs easily out to a redshift 2: ~ 10 
(|Blain et al.ll20o3 ), the fact that objects at 2 ^ 4.5 do not 
appear in our sample demonstrates that they do not exist 
in large quantities, and would therefore require much larger 
surveys to find them. Only if many of the identifications for 
these AzTEC sources are in fact more complicated (as in 
the case of HDF 850.1) may the door still be open for a sig- 
nificant fraction of the SMC population to lie at generally 
higher redshifts (2 > 4.5). 



4 SPECTRAL ENERGY DISTRIBUTIONS 

With redshift estimates in hand, we are now in a posi- 
tion to probe the rest-frame SEDs of our sample. Although 
we have photometry at a number of wavelengths spanning 
3.6 (i.m to 20 cm for most of the objects, the most inter- 
esting new constraints that we place on these SEDs is the 
shape of their rest-frame far-IR emission that peaks near 
100 \im in the rest-frame, produced by thermal dust grain 
emission. This emission accounts for most of the bolometric 
luminosity in SMGs, and it is generally believed to be pro- 
duce d by optically-ob scured star formation in most cases 
fe.g. lBlain eral]|2002l ). much like locally observed ULIRGs. 
The far-IR luminosity is therefore crucial for estimating star- 
formation rates. The far-IR SED also provides a direct probe 
of the total dust mass in a galaxy. However, both the bolo- 
metric luminosity and dust mass are critically dependent 
on the dust temperatur e, T, and the dust grain emissiv- 
ity, 13 jHildebrandl lioS^V Due to a dearth of data at the 
necessary wavelengths, spanning ~ 100-1000 ii.m, most au- 
thors either attempt to fit a simple 3-parameter modified 
blackbody spectrum for a population of dust grains at a 
single temperature. Si, = Av^ B^iT) (where A is the am- 
plitude), or adopt a single SED and normalize it to the 
(sub) mm data point. Since only a single (sub) mm data 
point is usually available, this latter compromise is often 
made. A census of recent studies finds broad agreement 
that the most typical values are Td = 30-35 K for SMGs, 
with an allowe d range t hat is som ewhat b roader than this 
JChapman et al. 2005; Kovacs et al. 20 06: Pope et al.|[2006l : 
iHuvnh et all I2OO7I : iCoppin et aLlbOOSI ). However, the esti- 
mates of Td and (3 are highly correlated, because of the lim- 
ited range of wavelengths for which data exist. 

In GOODS-N the combination of 1.1mm and 850 vim 
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flux densities sample wavelengths longward of the rest-frame 
far-IR peak. The ratio S'gso/S'i.i defines a family of 2- 
parameter SEDs (T and jS) for each source which we will 
use to check for consistency with previous measurements of 
the thermal SEDs of SMGs at typically shorter wavelengths. 



4.1 Correcting for flux density bias 

We first estimate un-biased 1.1mm and 850 y im flux den- 
sities for the AzTEC sources. As discussed in IPerera et al.l 
l|2008[ ). the 1.1mm flux densities are biased high because 
they are selected from a low SNR list of peaks coming from a 
counts distribution that falls steeply with increasing bright- 
ne ss. The correc tion for this bias followed the prescription 
of ICoppin et al. I (|2005i l. and we adopt those posterior flux 
density distributions here. While this correction does not ac- 
count for the additional effect of source blending, the AzTEC 
GOODS-N survey is shallower than the estimated confusion 
limit, and the sources that appear to be confused have been 
fit explicitly in this paper using two components. Rather 
than cross- matching the AzTEC catalogue with the SCUBA 
catalogue to obtain 850 ^m flux densities, which itself suffers 
flux density bias (and since the SCUBA data are also too 
shallow in some areas to provide flux densities for many of 
our sources), we instead directly measure the 850 ^.m map at 
the positions of proposed counterparts. Provided that these 
counterparts are correct, and 850 \im source confusion is neg- 
ligible, this photometry yields un-biased 850 (xm flux den- 
sities with symmetric Gaussian uncertainties for all 24/29 
AzTEC sources that land within the region of SCUBA cov- 
erage. However, due to the wide range in sensitivities only 
9 objects have 850 \im detections with a signiflcance of at 
least 3- a. 



4.2 Searching for '850 \x.m dropouts' 

It has been suggested that in regions where observations 
at both 850 Jim and ~l.lmm exist '850 (i.m dropouts', 
i.e. sources that are detected by AzTEC but not by SCUBA, 
can be used to select pr edominantly higher-redshift sources 
l|Greve et al.ll2004l2008l ). This technique is expected to work 
for the same reason that the AzTEC redshift distribution 
is slightly higher than the SCUBA sample: there is an in- 
creased submm negative Tf-correction at 1.1mm compared 
to 850 |i.m (i.e. the ratio of 850 jxm to 1.2 mm flux density, 
i5'85o/'S'i .2, decreases with redshift, seen for example in Fig- 
ure 4 of lEales et al.|[2003 ). In this study, we proceed by first 
testing the hypothesis of a single intrinsic observed fiux den- 
sity ratio R = Sgso/Si.!, by measuring R for several high- 
SNR objects selected in the AzTEC map, and then search- 
ing for dropouts in the SCUBA map relative to this average 
colour. We also repeat this analysis in the opposite direc- 
tion (for completeness), searching for SCUBA sources that 
are dropouts in the AzTEC map. 

We measure R using sources with deboosted fiux den- 
sities that have significances > 3-a in both bands: AzGN 1, 
AzGN 3, AzGN 7 and AzGN 8, giving similar values 1.72, 
1.96, 1.67 and 2.04 respectively. We adopt the mean R = 1.8. 
For reference, thermal emission from a galaxy with T = 30 K 
and /3 = 1.5atz = 2.5 would give an observed ratio 1.85. 
In order to compare our measurement with values reported 



for SCUBA (850 Jim) and MAMBO (1.2 mm) overlap, we 
use the same model SED to estimate how much larger the 
<S'85o/'S'i.2 ratio would be, finding that the scaled result is 2.3, 
near the cen tre of the distributions reported bv lGreve et al.l 
(j2004 [ 20081 ^ . Similarly we scale our result to estimate the ra- 
tio 5890/51.1 for 890 ^m SMA follow-up of AzTEC sources, 
finding a ra tio of 1.6. This value is consistent with 1.4 ± 0.3 
reported bv I Younger et all (|2007l '). 

Next we use our measured R to scale the deboosted 
1.1mm fiux density distributions for the entire sample to 
850 \im. For simplicity we approximate the scaled 850 (xm 
predicted fiux density distributions as Gaussians with mean 
values Sin given by the modes, and standard deviations am 
as half of the 68 per cent confidence intervals. If our hypoth- 
esis of a single observed ratio were true, given the observed 
850 \im data with mean fiux densities Sd and standard devia- 
tions (Td, we would expect the residuals (sm — Sd)/\/ o"m + o"| 
to be normally distributed with mean and standard devi- 
ation 1. For the 20 objects that do not have S-cr detections 
at SSOjim we calculate the sample mean and standard de- 
viation of the residuals, giving —0.1 ± 1.4. This calculation 
confirms that a ratio of 1.8 is a good estimate for the cen- 
tral value of the observed distribution Sgso/^i.i (left panel 
of Fig.[g. 

Since the residual is broader than expected (by a factor 
~ \/2), we conclude that the intrinsic spread in R produces 
uncertainties of the same order as our measurement errors. 
Note that the spread in 7? for a T = 30 is', /3 = 1.5 SED 
from redshifts z = 1-4 is only 2.1-1.6, so that part of the 
measured spread must be due to a range of rest-frame dust 
emission spectra in addition to the redshift distribution. We 
note that this scatter is roughly symmetric: the 850 \ira mea- 
surements fall above the expected values about as often as 
they fall below (in the left panel of Fig. |6} . 

These calculations are biased to 1.1mm selected 
sources, and could in principle be different for an 850 \ira 
selected catalogue. We therefore repeat the procedure, start- 
ing with deboosted fiux densities for the 20 robustly iden- 
tified SCUBA sources mentioned in Section [31 scaling them 
to 1.1 mm by dividing them by R, and comparing these pre- 
dictions to photometry at the locations of their counter- 
parts in the AzTEC map (right panel of Fig. In this 
case the scatter is clearly asymmetric, with a number of 
the SCUBA sources appearing fainter than expected in the 
AzTEC map. However, these outliers are primarily low-SNR 
SCUBA detections (indicated by the large horizontal error 
bars), in particular GN 3, GN 5, GN 7, GN 16 and GN 22 
which have uncert ainties rang i ng fro m 1.5-4.5 mjy. We hy- 
pothesize that the lPope et al. (|2006l ) deboosting recipe did 
not sufficiently correct these sources, i.e. they should have 
been shifted further to the left in this plot. This explanation 
is plausible because their deboosting factors were extrapo- 
lated from those cal culated for the SCUBA SHADES survey 
ICoppin et al.l l2006). which were derived for sources with a 
different noise distribution. Furthermore, the fact that the 
AzTEC selected source list does not exhibit this problem 
(left panel of Fig. |6} also points to an issue with the SCUBA 
deboosting calculation, rather than the map itself. 

Finally, we use the AzTEC and SCUBA selected source 
lists to plot histograms of the measured 585o/5i.i flux ratios 
m Fig. [3 Here we have convolved each source with its uncer- 
tainty distribution before adding it to the total histogram. 
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Figure 6. Measured flux densities (vertical axes) compared with predicted flux densities (horizontal axes) for AzTEC (1.1mm) and 
SCUBA (850 i^m) selected source lists, left: Measured SCUBA flux densities at the locations of counterparts proposed for AzTEC sources. 
The predicted SCUBA flux densities were scaled from the deboosted AzTEC measurements using the ratio of iSgso/Si.i fit to the 4 
most significant simultaneous detections in both the AzTEC and SCUBA maps (thick square symbols), yielding 1.8. We also plot the 
relation for the remaining 20 objects measured in both bands (diamonds). The symbols and error bars are coloured more lightly if they 
are less significant (combined observed and predicted 850 |^m uncertainties). This plot shows symmetric scatter about the single-colour 
relation, and hence no clear evidence for 850 (xm dropouts that would fall systematically below the line, right: AzTEC fiux densities at 
th e locations of counterparts to SCUBA sources. The predicted values are derived from the deboosted SCUBA fiux densities reported 
in lPope et al.l |200G) by dividing them by the same factor of 1.8. The low outliers are primarily low-SNR SCUBA measurements (larger 
horizontal error bars) suggesting that they have been insufficiently deboosted. 



SO that higher-SNR measurements are more sharply peaked 
and therefore contribute more to the shape. In the left panel 
of this figure the resulting histograms for each sample are 
shown to be nearly indistinguishable, with modes that are 
coincident with the value R = 1.8 measured for the highest- 
SNR detections in the AzTEC map (vertical dotted line). 
Since we also have redshift information, we split the sam- 
ples into objects above and below z = 2.5, allowing us to 
search for a systematic trend. For the AzTEC selected sam- 
ple, there is no significant difference in the ratios for low 
and high-redshift objects. On the other hand, the mode of 
the SCUBA sample shows a mild trend to higher values of 
R with increasing redshift. However, given the evidence for 
incorrect deboosting of the SCUBA flux densities mentioned 
above we do not believe this trend is significant. 

From these measurements we therefore conclude that 
there is no evidence for 850 [im dropouts in the 1.1 mm map. 
In fact, the observed trend in the Ssso/Si.i colour tends 
to go in the opposite direction: AzTEC fiux densities ap- 
pear fainter than they should be starting with a SCUBA 
selected list (but not vice versa). How then can this result 
be reconciled with the measured shift in the AzTEC redshift 
distribution compared to SCUBA (Fig. ^ , which can only 
be a result of enhanced brightness at 1.1mm at larger val- 
ues of z? Given the low-SNR measurements used here, the 
wavelength proximity of the 850 vim and 1.1 mm filters, and 



the fact that t here appear s to be a bias in the SCUBA flux 
densities from lPope et al.l (|2006l ). we believe that trends in 
the observed colour are completely obscured by uncertain- 
ties. Much deeper (and higher spectral resolution) observa- 
tions would be needed to detect differences in the colours of 
850 [im and 1.1 mm source populations. 



4.3 Constraints on T and l3 

Given the T-/3 degeneracy, low SNR, and wavelength prox- 
imity of the SCUBA 850 ^m and AzTEC 1.1 mm filters, the 
SED constraints for individual objects are extremely noisy. 
However, with this data-set it is possible to examine the con- 
straints that we can place on the family of T and /3 consistent 
with the sample. Note that this is different to fitting SED 
templates based on objects in the local Universe using all 
of the available radio -IR photometry (e.g. iPope et"aL I I2OO6I : 
iMagnelh et"aLll2009^ . The purpose of such flts is to give the 
best estimates of source properties, such as bolometric lu- 
minosities and star-formation/ AGN fraction. However, the 
goal in this paper is to provide un-biased measurements that 
in the future may be used to constrain redshift-dependent 
SED libraries. 

We have extracted a subset of the galaxies for which 
there are robust counterparts, 850 ^.m flux densities with a 
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Figure 7. Histograms of the R = Sg5o/5i,i flux ratio. In each case the vertical dotted hne corresponds to a value of R= 1.8 calculated 
from high-SNR, sources in Fig |6] The top horizontal axes assume a /3 = 1.5, T = 30K modified blackbody spectrum to predict the 
corresponding Ssso/'Si.2 colours for comparison with MAMBO results, left: The black and grey histograms show the distributions of 
AzTEC and SCUBA selected sources, respectively. As in Fig. [5] the uncertainty of each object is convolved with its error distribution 
before placing it in the histogram, such that higher-SNR measurements are more sharply peaked, right: Here the samples have been 
divided into subsets above and below redshift z = 2.5. Although similar numbers of objects lie in each redshift bin, the higher redshift 
bin is dominated by the high-S/N detections of AzGN 1 and AzGN 1.2. 



significance of at least 3-cr, and spectroscopically measured 
redshifts: AzGN 1, AzGN 1.2, AzGN 7, and AzGN 16. We 
then plot the joint likelihood surface of the observed 850 vim 
and 1.1mm flux densities given a range of models with T 
and 13 common to each object (marginalizing over the am- 
plitudes A for each source as they are irrelevant to the joint 
distribution) in Fig. [8] This surface, as expected, produces a 
long anti-correlated valley between T and /3. For comparison 
we also plot the locus of SED models that would produce 
the sa me peak S{v) as the best-fit model from lCoppin et ahl 
l|2008l l, T = 31K assuming 13 = 1.5. As they use shorter- 
wavelength 350 |J.m follow-up of SMGs, the constraints are 
somewhat orthogonal to those described in this paper, show- 
ing the ability of observations at those wavelengths to reject 
the highest-temperature and lowest-/? SEDs. 

There is some evidence that this AzTEC sample lies 
at slightly higher redshift (Fig. [SJ , but the fact that we do 
not find significant '850 \i.m dropouts' suggests that, within 
the precision of the data, the SEDs are not different. At 
face value. Fig. [S] suggests that the best-fit value at the in- 
tersection of the valley in the contour plot with the solid 
line occurs at T ~ 30K, and (5 ~ 1.75. However, it is clear 
that a large family of SEDs with temperatures ranging from 
T ~25-40K and corresponding (3 ~2-l are easily allowed. 



5 CONCLUSIONS 

We have used the rich multi-wavelength data-set in GOODS- 
N to identify robust radio and IR counterparts for sources 
detected at 1.1mm using AzTEC. Of the 29 sources, only 
1 or 2 are expected to be false positives. We find robust 
counterparts for 22 objects (false identification probabilities 
P < 0.05), although one object {AzGN U/HDF 850.1) is 
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Figure 8. Dotted lines are contours in the joint likelihood func- 
tion (0.6, 0.1 and 0.001 times the maximum likelihood) demon- 
strating the family of values for /3 and T consistent with a si- 
multaneous fit to the 5 galaxies for which spectroscopic redshifts 
have been obtained, and there exist at least 3-a measurements 
of the 850 \xm fiux densities: AzGN 1, AzGN 1.2, AzGN 3, 
AzGN 7, and AzGN 16. The simple 3-parameter rest-frame SED 
Aul^ Bi,{T) is redshifted, and compared to the observed 850 \Jtm 
and 1 . 1 mm fiux densities for each object with a grid of parameter 
values, and then marginalized over A. The final two-dimensional 
likelihood surface is the product of the five individual marginal- 
ized distributions. The solid line is the locus of SED s that produce 
the same peak in S(i') as the best-fit model from lCoppin et al.l 
(200l), T = 31K and /3 = 1.5 such that we can see the effect 
of an additional shorter-wavelength measurement to reject the 
highest-temperature / lowest-/? fits. 
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dropped from the sample due to confusion about its identifi- 
cation, and tentative associations for the remaining 8 objects 
are also provided. These counterparts have an astrometric 
precision of ~ 1 arcsec, a significant improvement given the 
18 arcsec FWHM AzTEC beam and low SNR. 

We find spectroscopic redshifts for 7 of the robustly 
identified sources in the literature, and provide photomet- 
ric redshifts or limits for the remaining objects. Restrict- 
ing ourselves to the 18 objects with robust counterparts 
and redshift estimates (spectroscopic or photometric, ex- 
cluding limits), we measure a median z — 2.7, with an in- 
terquartile range 2.1-3.4. The 850 [im sources in this field, 
selected in a similar way, have a median redshift z — 2.0 
with an interquartile range 1.3-2.6. We use K-S and M-W 
non-parametric tests to evaluate the significance of this shift 
to higher redshift in the 1.1mm map, finding chance prob- 
abilities pks = 0.05 and pmw ~ 0.01 that both surveys 
sample the same redshift population. Given the large uncer- 
tainties in individual photometric redshifts, we used Monte 
Carlo simulations to evaluate the spread in probabilities pro- 
duced by the two tests consistent with our samples, finding 
Pks < 0.15 and pmw < 0.04 at a confidence level of 68 per 
cent. 

For the entire overlapping region between the SCUBA 
and AzTEC maps, we perform un-biased flux density mea- 
surements in the SCUBA map at locations of identiflcations 
for AzTEC sources. Using the 4 most significant (3-a) de- 
tections in the two maps we find a mean observed 850 ^im 
to 1.1mm flux density ratio Sgso/Si.! = 1.8. For the re- 
maining 20 sources, we observe a symmetric scatter in the 
observed ratio which appears to be produced in equal quan- 
tities by intrinsic spread in spectral properties, and measure- 
ment noise. We also examine the ratios Ssso/Si.i for objects 
selected at 850 \im, finding that they are also generally con- 
sistent with this value, although it appears that the lower- 
SNR 850 \im flux densities may be biased high. Finally, we 
unsuccessfully searched for trends in this flux density ratio 
with redshift for both samples. We therefore do not see evi- 
dence for 850 ^.m dropou ts in the 1.1 mm map as reported in 
ICreve et al"] l|2004l2008l ). While we believe that such a trend 
must exist in the underlying SMC population to produce 
the mild differences in the redshift distributions mentioned 
above, it is undetectable when comparing ~4-cr detections 
in wavelength bands that are so close to eachother. 

We test the hypothesis of a single temperature, T, and 
dust emissivity index, (3, for the ensemble of sources having 
robust identifications and photometric redshift estimates. 
Given the degeneracy between these parameters (since we 
have only two photometric measurements at different wave- 
lengths), we assume the same mean rest-frame far-IR peak 
as found in other studies, finding that T = 30 K and /3 — 1.75 
are consistent with all of the data. However, given the SNR, 
these measurements still provide only a weak constraint, and 
data at shorter rest-frame far-IR wavelengths would be re- 
quired to tighten up the allowable range of SEDs. SCUBA-2 
450 ^lm, as well as S PIRE and new BLAST 250, 350 and 
500 ^m surveys (e.g. iDevhn et all l2009l : iDve et all |2009| ) 
should be particularly useful. 



Table Al. New source positions and raw 1.1 mm map flux densi- 
ties resulting from simultaneous two-source fits. Source Az 1.2 is a 
new object in this paper, w hereas the other thre e were originally 
detected in AzTEC maps in lPerera et al.l l|2008l ). 



AzTEC 


R.A. 


Dec. 


Map Flux Density 


ID 


m s "J 


(° ' ") 


(mjy) 


1 


12 37 11.99 


+62 22 11.1 


12.73 + 0.99 


1.2 


12 37 09.15 


+62 22 02.1 


4.14 + 0.98 


14 


12 36 52.23 


+62 12 25.2 


4.11 + 0.97 


22 


12 36 49.12 


+62 12 13.1 


4.56 + 0.97 
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APPENDIX A: DATA TABLES 

Here we provide proposed identifications and multi- 
wavelength photometry for all of the sources. Table lAll 
provides updated positions and raw map flux densities for 
AzTEC 1.1 mm sources that required deblending. Tables 1X21 
and lA3l summarize the identifications and SEDs respectively. 



APPENDIX B: NOTES ON EACH SOURCE 

This section gives detailed information on each source not 
provided in the tables of Appendix 1X1 

Bl Robust Identifications 

Sources in this category have potential counterparts with 
P < 0.05 within 6 arcsec. 

AzGN 1: GN 20 from jPope et all (|2006| ). This source 
has been deblended from AzGN 1.2 (see Table IA1|) . and 
the submm emis sion was also localized using the SMA 
(|lono et al.l |2006| ). A spectr oscopic reds h ift of 4.055 for 
this source was reported in iDaddi et al.l (|2008l ) based on 
molecular CO emission detected with the IRAM PdBL, 
with some confirmation based on optical spectroscopy in 
Pope et. al (submitted). 

AzGN 1.2: This was already known to be a second compo- 
nent of AzGN 1 from the SCUBA data. In the AzTEC map 
the source was deblended by performing a simultaneous fit 
of two scaled effective PSFs using the peak at the location 
of Az GN 1 and the position of GN 20.2 from jPope et all 
(|2006l ) as starting values (see Table IXT]) . Positions and flux 
densities were then allowed to vary (a total of 6 parame- 
ters). The reduced for this fit decreased to 0.99 from 1.15 
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Table A2. Radio and Spitzer identifications of AzTEC sources (procedure described in Section [2)l. Counterpart distances in brackets 
employed a 10 arcsec search radius. P values in boldface emphasize robust counterparts with values < 0.05. Spectroscopic redshifts are 
given in the column labelled 2apoc (references for these measureme nts given in Appen dix IB )| . Photometric redshifts based on Spitzer IR 
flux densities from Table {K3\ are calculated using Equation 2 from lPope et aU ||2006| ) and given in the penultimate column, 2;^. These 
redshifts have uncertainties Azj^ = 0.15(1 + 2), and are biased low at z > 3. Photometric redshifts based on the (sub)mm-to-radio colours 
are given in the last column, Zrm- The quoted 68 per cent uncertainties are theoretically derived; in this paper we assume an empirically 
measured symmetric error Azmi = 0.8. 



AzTEC SCUBA Radio Spitzer Redshift 



ID 


ID 




R.A. 


Dec. 


Dist. 




R.A. 


Dec. 


Dist. 


P 






Zrm 








m s J 


{° 




") 


(") 




m 


(° 




") 


{") 








1 


20 


12 


37 


11.88 


+62 


22 


11.8 


1.0 


12 


37 


11.88 


+62 


22 


12.1 


1.3 


0.003 


4.055 


2.7 


3.»_Q 7 


1 o 


on o 


12 


37 


08.78 


+62 


22 


01.8 


z.D 


1 o 


37 


08.77 


+62 


22 


01.8 


9 7 
Z. / 


U.UUo 




9 ^ 
Z.O 


q 1+1.2 

•-'■^-0.2 






12 


37 


09.73 


+62 


22 


02.5 


4.1 


12 


37 


09.57 


+62 


22 


02.1 


2.9 


0.037 




3.1 


9 t-+0.8 
^■^-0.8 


2 




12 


36 


31.93 




17 


14.7 


5.3 


12 


36 


31.92 




1 7 
J- 1 


14 6 


5.2 


0.051 




3.2 


2.4j;2-3 


3 


10 


12 


36 


33.42 


+62 


14 


08.7 


0.6 


12 


36 


33.40 


+62 


14 


08.4 


0.6 


0.002 


4.042 


2.3 


n -1+1.4 

■^■^-0.2 


4 




12 


35 


50.26 


+62 


10 


41.3 


3.1 


12 


35 


50.35 


+62 


10 


41.8 


2.7 


0.030" 




2.9 


2 4+2-2 


5 




12 


37 


30.78 


+62 


12 


58.7 


2.6 


12 


37 


30.75 


+62 


12 


58.4 


2.2 


0.007 




2.0 




6 




12 


36 


27.26 


+62 


06 


05.7 


1.6 


12 


36 


27.21 


+62 


06 


05.7 


1.2 


0.007 




3.0 




7 


39 


12 


37 


11.32 


+62 


13 


30.9 


4.4 


12 


37 


11.34 


+62 


13 


31.0 


4.3 


0.014 


1.996 


1.7 


9 7+0.9 
^■'-0.8 






12 


37 


11.99 


+62 


13 


25.6 


4.5 


12 


37 


11.99 


+62 


13 


25.7 


4.4 


0.032 


1.992 


2.0 


>2.3 


8 


12 


12 


36 


46.04 


+62 


14 


48.6 


(6.9) 


12 


36 


46.07 


+62 


14 


48.8 


(7.0) 


0.037 




2.0 


3.0in 






12 


36 


46.80 


+62 


14 


45.3 


(7.5) 


12 


36 


46.88 


+62 


14 


47.2 


(9.0) 


0.050 


2.006 


0.7 




9 


37 


12 


37 


38.16 


+62 


17 


37.0 


1.6 


12 


o 1 


oo.zu 


+62 


17 


36.4 


0.9 


0.013 


3.1900 


2.4 


>3.0 


10 




12 


36 


27.54 


+62 


12 


17.8 


3.5 


12 


OD 


97 /I fi 


+62 


12 


18.0 


3.1 


0.066" 
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9 1+2.2 
7 


11 


27 


12 


36 


35.89 


+62 


07 


03.8 


3.1 
















0.027 






2.8+}'i 

—0.2 


12 




12 


36 


32.65 


+62 


06 


21.1 


4.7 


12 


36 


32.65 


+62 


06 


21.3 


4.9 


0.047 




2.8 


>1.3 


13 




12 


35 


54.23 


+62 


13 


43.8 


2.9 


12 


35 


54.28 


+62 


13 


43.4 


3.3 


0.033" 




2.2 


2 6+0.9C 


14 


14 


12 36 52.07'' 


+62 


12 25. 7^* 


1.2 






















15 




12 


35 


47.93 


+62 


15 


29.2 


5.0 


12 
12 


35 
35 


48.09 
47.87 


+62 
+62 


15 
15 


29.3 
28.2 


3.9 
5.7 


0.018 






>1.6 


16 


04 


12 


36 


16.09 


+62 


15 


13.8 


4.3 


12 


36 


16.10 


+62 


15 


13.6 


4.5 


0.038 


2.578 


2.0 


2 4+2-1 






12 


36 


15.80 


+62 


15 


15.1 


4.0 


12 


36 


15.82 


+62 


15 


15.4 


3.6 


0.039 




4.5 


n Q+0.3 


17 




12 


35 


39.92 


+62 


14 


42.1 


(7.6) 


12 


35 


39.95 


+62 


14 


40.8 


(6.5) 


0.062 






>1.7 






12 


35 


39.91 


+62 


14 


30.8 


(7.1) 


12 
12 


35 
35 


39.97 
39.76 


+62 
+62 


14 
14 


30.7 
30.7 


(6.9) 
(7.9) 


0.067 






>1.7 
>1.7 


18 


38 


12 


37 


41.16 


+62 


12 


20.5 


3.8 


12 


37 


41.16 


+62 


12 


20.9 


3.5 


0.032 




2.1 


2 4+0.6 
^■^-2.0 






12 


37 


41.63 


+62 


12 


23.6 


5.8 


12 


37 


41.66 


+62 


12 


23.6 


6.0 


0.051 




2.0 


>2.2'= 


19 




12 


36 


04.40 


+62 


07 


02.7 


2.6 


12 


36 


04.43 


+62 


07 


02.9 


2.8 


0.022 






>1.5 


20 




12 


37 


12.48 


+62 


10 


35.4 


3.0 


12 


37 


12.51 


+62 


10 


35.6 


2.8 


0.030 




3.3 


>2.3 


21 




12 


38 


00.80 


+62 


16 


11.7 


1.5 


12 


38 


00.81 


+62 


16 


11.7 


1.4 


0.016" 




2.3 


>1.6 


22 




12 


36 


48.60 


+62 


12 


16.1 


4.7 


12 


36 


48.65 


+62 


12 


15.7 


4.2 


0.083" 




1.7 


>0.7 


23 




12 


37 


16.67 


+62 


17 


33.2 


1.4 


12 


37 


16.67 


+62 


17 


33.3 


1.4 


0.001 


1.1460 


1.6 


>1.8 


24 


23 


12 


36 


08.58 


+62 


14 


35.3 


(6.4) 


12 


36 


08.60 


+62 


14 


35.3 


(6.4) 


0.079 




2.7 


2 4+2-0 


25 




12 


36 


51.72 


+62 


05 


03.0 


4.1 
















0.021 






>0.7 


26 


40 


12 


37 


13.86 


+62 


18 


26.2 


0.6 


12 


37 


13.85 


+62 


18 


26.2 


0.5 


0.000 




2.6 


>0.7 


27 




12 


37 


19.62 


+62 


12 


20.3 


1.3 


12 


37 


19.61 


+62 


12 


20.9 


0.9 


0.034* 




3.1 


>0.7 






12 


37 


20.01 


+62 


12 


22.0 


2.1 


12 


37 


19.99 


+62 


12 


22.6 


2.2 


0.050* 


2.4600 


2.0 


>0.7 


28 




12 


36 


44.03 


+62 


19 


38.8 


4.2 


12 


36 


44.03 


+62 


19 


38.4 


3.9 


0.071" 




2.5 


>0.7 



"Identified in 3-(t radio catalogue 
^Identified in MIPS 24 yxm catalogue 
'^Photometric redshift calcula tion excludes e xtremely noisy 850 ^m photo metry 
''Position for HDF 850.1 from lDunlop et al. , (2004). See also Co wie et al.l (|2009). 



when only a single source was used, justifying the addition 
of the extra parameters. This procedure yields a clear 4.2- 
a sou rce which corresponds to GN 20.2 from IPope et al.l 
l|2006l ). There are two possible radio identifications, one 
fainter object is 4.1 arcsec away with P — 0.037, and the 
other brighter source is 2.6 arcsec away with P = 0.005. The 
brighter object is t he claimed counterpart to AzGN 1.2 from 
iPope et al.l l|2006l ). However, both p otential coun t erpart s 
mentioned here are also discussed in iDaddi et al.l l|2008h . 



They detect a faint emission line from the brighter radio 
source which they argue to be consistent with molecular CO 
emissi on at a redshift 4.051. We concur with iDaddi et al.l 
l|2008t ) that the brighter object is likely an AGN, based on 
its relatively large radio/mm flux ratio. 

AzGN 3: The AzTEC detection, and proposed coun- 
terpart, are both coincident with the S CUBA source 
GN 10 and identification f rom iPope et al.l (|2006| ). Similar 
to AzGN 1 and AzGN 1.2. lDaddi et al.l l|2009l )ldentified a 
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Table A3. Photometry data list ed in the same o rder as the identifications in Table IA2I The 850 \xm flux densities are measured 
directly from the SCUBA map of IPope et al at the positions of each proposed counterpart, and should therefore be un-biased 

at each correct pos ition. 24 [om upper-limits are given at a significance of 3-cr. De-boosted AzTEC 1.1mm flux densities are taken from 
iPerera et al] |j200^) except for AzGN 1, AzGN 1.2, AzGN 14 and AzGN 22, which are calculated using the raw map flux densities from 
Table lAl] but adopting the same prior. 



AzTEC 


1.4 GHz 


1.1 mm 


850 |am 


24 |am 


8 |am 


5.8 ^mi 


4.5 ^mi 


3.6 |am 


ID 


(^Jy) 


(mjy) 


(mJy) 




(MJy) 


(/^Jy) 


(^Jy) 


(^Jy) 


1 


72.8 ±12.3 


11.81±+;n? 

— l.U ( 


20.30 ±2.10 


65.5 ±4.5 


25.3 ± 1 


.5 


15.9 ± 1, 


.3 


9.2 ±0.8 


6.8 ±0.6 


1.2 


173.5 ±6.3 


2.87±+V?= 

— 1 .ZD 


9.90 ± 2.30 


20.2 ± 3.5 


9.4 ± 1 


.1 


6.1 ± 1, 


,0 


3.9 ±0.5 


3.9 ±0.6 




35.2 ±6.3 


9.90 ±2.30 


12.9 ± 3.1 


14.9 ± 1 


.5 


11.8 ± 1, 


,3 


9.8 ±0.8 


8.0 ±0.9 


2 


26.2 ±4.6 


5 9i±+i-02 


5.74 ±2.81 


37.7 ±4.8 


24.9 ± 1 


.5 


17.0 ± 1, 


.3 


10.9 ±0.8 


6.4 ±0.6 


3 


36.0 ±4.2 


5.35±+?-g^ 


10.50 ± 1.59 


26.2 ± 5.1 


5.1 ± 1, 


.1 


2.7 ±0, 


.9 


2.0 ±0.4 


1.2 ±0.4 


4 


34.1 ± 10.7 


4.69±+l°5 
— i.ud 




92.5 ±4.2 


37.0 ±1, 


.7 


26.0 ± 1, 


.5 


16.0 ±0.8 


9.6 ±0.9 


5 


128.2 ± 8.0 


4.13±+3,°^ 

— 0.98 


4.37 ±3.92 


181.0 ± 6.5 


28.9 ±1, 


.7 


18.9 ± 1, 


,3 


12.7 ±0.8 


9.2 ± 0.9 


6 


58.5 ± 13.1 


^^ — 1.00 




11.3 ± 8.3 


9.9 ±1, 


.1 


4.1 ± 1, 


,0 
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spectroscopic redshift of 4.042 using CO emission detected 
with the IRAM PdBI. 

AzGN 4: There is a single unambiguous radio counterpart. 

AzGN 5: There is a single unam biguous rad i o cou nter- 
part. This object is not part of the lPope et al.l (|2006l ) cat- 
alogue (probably because of the high noise in this region), 
but was originally discovered in the SCUBA jiggle map of 
IWang et all (|2004l ') (GOODS 850-6). Our proposed counter- 
part corresponds to object 'c' from their analysis. We note 
a large discrepancy between the 850 [im flux density of this 
source at the coordinates of the proposed counterpart in 
the HDF supermap, 4.4 ± 3.9 mJy, used in thi s paper (Ta- 
ble IA3II . compared with 19.41 ± 3.2 mJy from IWang et all 
(I2OOJ) which is due at least in part to the lack of correction 



for Eddington bias in that work. Based on the 1.1mm flux 
density and our measured ratio Ss5o/Si.i — 1.8, we would 
expect 5850 — 7mJy. 

AzGN 6: There is a single unambiguous radio counterpart. 

AzGN 7: This obje ct is also a s ignific ant SCUBA source, 
GN 39, described in IWaU et all (l2008t l. We find the same 
two radio identificatio ns with spe ctroscopic redshifts 1 .996 
(|Chapman et al.ll2005l ) and 1.992 (|Swinbank et al.ll2004^ re- 
spectively. 

AzGN 9: This objects i s also a significant SCUBA source, 
GN 37 from iPope et al.l (|2006l ') for which we fi nd the same 
ident ification with spectroscopic redshift 3.190 (jCowie et al.l 
l2004h . 

AzGN 11: This object is close to the SCUBA source 
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GN 27 from IPope et all (120061 ) that is classified as 'ten- 
tative' since it does not exceed their threshold for posterior 
deboosted flux density probability above 0. There appear to 
be three distinct peaks in the radio map that fall between the 
SCUBA and AzTEC peaks (but closer to the AzTEC posi- 
tion). Two of these peaks are coincident with Spitzer sources 
at redshifts 0.56 ar id 0.276. The ACS /IRAC counterpart to 
GN 27 claimed in IPope et al.l (120061 ) is to the south of the 
three objects that we find, and is not a significant source in 
our 1.4 GHz map. The third source is considerably fainter 
and redder in the IRAC channels. However, due to blending 
with the other sources we are unable to produce the SED 
at Spitzer wavelengths in Table IA3I We rule out the low- 
z identifications based on their highly unlikely SEDs which 
would exhibit extraordinarily low radio/(sub)mm flux den- 
sity ratios. 

AzGN 12: There is a single unambiguous radio counter- 
part. 

AzGN 13: There is a single unambiguous radio counter- 
part within 6 arcsec. We note that there are two radio 
sources 7.9 arcsec (12:35:53.90 -1-62:13:37.1) and 8.7 arcsec 
(12:35:53.24 -1-62:13:37.5) to the south that are coincident 
with a low-significance peak o f SCUBA emission. Th e sec- 
ond source was also part of the lChapman et al.l (|2005h sam- 
ple. While they are not considered likely in our analysis, we 
list them here for completeness. Had we used a search radius 
of 10 arcsec they would hav e P values of 0.103 and 0.108, 
and he at redshifts . 8770 (iReddv et all |2006| ) and 2.098 
l|Chapman et al.ll2005l : [Reddv et al.ll2006f ). respectively. 

AzGN 14: This object is H DF 850.1, the highe st SNR 
SMG from the HDF-N map of HugheseTaD l|l998h . Simi- 
lar to the case for AzGN 1 and AzGN 1.2, this source has 
been de-blended from AzGN 22 (see Table IA1|) . Various 
counterparts have been suggest ed, with two det ailed stud- 
ies giving different answers. D unlop et all (|2004 ) subtract 
the emission from an elliptical galaxy that is believed to ob- 
scure a faint background K-hand counterpart at 12:36:52.07 
-1-62:12:25.7 (this p osition is given in Table IA2[I . Recently, 
ICowie et al.l (I2009D used the SMA to localize the source of 
the submm emission to 12:36:51.99 -1-62:12:2 5.83, conclud- 
ing th at this position is incompatible with the lDunlop et all 
l|2004h counterpart. For the purposes of measuring un-biased 
flux densities in the SCUBA and AzTEC maps, either po- 
sition may be used as the differences are small compared to 
the (sub)mm beams. Had we applied our counterpart search 
blindly we would have instead chosen a more distant radio 
source at 12:36:51.72 -1-62:12:21.36 (P = 0.036), although we 
would have rejected this counterpart based on the unrealistic 
SED. This case demonstrates the potential for counterparts 
identified with low values of P to be misleading. No Spitzer 
fiux densities are quoted for this special source, and it is 
excluded from most of the analysis in this paper. 

AzGN 15: There is a single unambiguous radio counter- 
part, although it appears to resolve into two objects in the 
Spitzer catalogue. Inspection of the map shows, however, 
that these data lie at the extreme edge of the Spitzer cover- 
age, and we suspect that the two sets of photometry are in 
fact for the same source. 

AzGN 16: Similar to A zGN 7, this objec t is also detected 
with SCUBA {GN 4 from lPope et al.ll2006l '). and has two ra- 
dio counterparts. The mor e distant counterpart h as a spec- 
troscopic redshift of 2.578 (| Chapman et al.ll2005l ). 



AzGN 18: This object i s also a significant SCUBA source, 
GN 38 from iPope et al.l l|2006h . and we find the same two 
radio counterparts. One of these counterparts technically 
misses the significance threshold with P — 0.051, but lies 
very close to the centre of the SCUBA emission peak. Con- 
versely, the other object, with P = 0.032, lies closer to the 
peak of the AzTEC emission. The photometric redshifts for 
both sources are consistent with being at the same distance. 

AzGN 19: There is a single unambiguous radio counter- 
part, and the position is coincident with a region of faint 
SCUBA emission. It falls off the edge of the MIPS 24 |i.m cov- 
erage, but within the IRAC footprint. There is a bright op- 
tical source within 1 arcsec of the radio position at z = 0.65. 
However, we reject this redshift and SED combination as it 
would imply an extraordinarily low radio/(sub)mm fiux ra- 
tio. Therefore the radio identification proposed here may be 
incorrect and AzGN 19 is associated with the \ow-z optical 
source instead, or it is correct and unrelated to the optical 
source. Since the (sub)mm/radio photometric redshift esti- 
mate is consistent with a more typical SMG (zrm > 1.5) we 
prefer the first hypothesis. 

AzGN 20: There is a single unambiguous radio counter- 
part, and the position is coincident with a region of faint 
SCUBA emission. 

AzGN 21: There is a single unambiguous radio counter- 
part. 

AzGN 23: There is a single unambiguous radio counter- 
part. This object has the lowest spectr oscopic redshift in 
the sample, z = 1.146 (|Cowie et al.ll2004l '). 

AzGN 25: This object has a single clear radio counterpart, 
but lies beyond the Spitzer IRAC coverage. The MIPS 24 \Lm 
image exhibits emission at the location of the radio counter- 
part, but is not part of the Spitzer catalogue produced by 
the GOODS team. 

AzGN 26: This object is also a significant SCUBA source, 
GN 40 from lWaU et al.l (120081 ). There is a single unambigu- 
ous radio identification. 

AzGN 27: This source is coincident with faint SCUBA 
emission. There are no radio sources with P < 0.05, even 
searching out to 10 arcsec. However, there are two MIPS 
24 (i.m sources 0.9 arcsec and 2.2 arcsec away, with P values 
0.034 and 0.05, respectively. The more distant source has a 
spectroscopic redshift z = 2.460 (Stern et al. in preparation), 
and photometric redshift estimates for the other source are 
broadly consistent with this value. As this object appears 
to have an extraordinarily faint radio flux density given the 
24 (i.m measurement, we warn the reader that these identifi- 
cations may be suspect. 



B2 Other Identifications 

For the remaining sources we give our best estimates for 
the counterparts, considering objects with 0.05 < P < 0.10 
and/or searching out to 10 arcsec from the AzTEC positions. 

AzGN 2: There is a radio source 5.3 arcsec to the north 
that slightly misses the significance cut (P — 0.051). This 
proposed counterpart is in the direction of a faint 2-a signal 
in the SCUBA map. 

AzGN 8: This object is also a significant SCUBA source, 
GN 12. There are two potential radio counterparts 7 arcsec 
and 9 arcsec away. Had we used a search radius of 10 arcsec, 
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the values of P would have been 0.037 and 0.050. Since the 
closer object would then have a lower P, and is also closer to 
the SCUBA peak, we prefer it as the most likely candidate, 
although strictly based on the 1.1mm and 1.4 GHz emis - 
sion they are both reasonable candidates. iPope et al.l (|2006l ) 
were unable to obtain an optical spectrum for the nearer 
candidate. The other candidate was detected serendipitously 
with Keck at a redshift z = 2.006. This second object was 
rejected by iPope et"aL I (|2006 l based on its distance from 
the SCUBA centroid, and the unlikely rest-frame UV emis- 
sion that would be uncharacteristic of a dusty star-forming 
galaxy. 

AzGN 10: This object has a single potential radio coun- 
terpart within 6 arcsec, and is also coincident with faint 
SCUBA emission. 

AzGN 17: This object is at the edge of the MIPS 24 ^m 
coverage, but within the IRAC footprint. There are two po- 
tential radio counterparts at distances 7.1 arcsec and 7.6 arc- 
sec. 

AzGN 22: This source lies on a region of extended emis- 
sion to the south-east of AzGN 14, which is also seen in 
the SCUBA map. To reduce the effect of blending with 
AzGN 14 we have performed a simultaneous 2-object fit, as 
with AzGN 1 and AzGN 1.2 (resulting positions and map 
flux densities given in Table IXT)) . The most likely counter- 
part, given the resulting centroid, is a radio source 4.7 arcsec 
away. 

AzGN 24: There is a single potential radio counterpart 
6.4 arcsec to the south. We consider this a likely identifi- 
cation because it falls be tween the AzTEC peak, and the 
SCUBA source GN 23 (|Pope et al.ll2006h . for which the 
same identification was proposed. 

AzGN 28: This object has a single potential radio coun- 
terpart 4.2 arcsec away, and is coincident with faint SCUBA 
emission. It has a redshift estimate zir — 2.5 
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